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ABSTRACT 



Aims. By focusing on the oscillations of the cross-sectional area and the intensity of magnetic waveguides located in the lower solar 
atmosphere, we aim to detect and identify magnetohydrodynamic (MHD) sausage waves. 

Methods. Capturing several series of high-resolution images of pores and sunspots and employing wavelet analysis in conjunction 
with empirical mode decomposition (EMD) makes the MHD wave analysis possible. For this paper, two sunspots and one pore (with 
a light bridge) were chosen as representative examples of MHD waveguides in the lower solar atmosphere. 

Results. The sunspots and pore display a range of periods from 4 to 65 minutes. The sunspots support longer periods than the pore - 
generally enabling a doubling or quadrupling of the maximum pore oscillatory period. All of these structures display area oscillations 
indicative of MHD sausage modes and in-phase behaviour between the area and intensity, presenting mounting evidence for the pres- 
ence of the slow sausage mode within these waveguides. 

Conclusions. The presence of fast and slow MHD sausage waves has been detected in three different magnetic waveguides in the 
lower solar photosphere. Furthermore, these oscillations are potentially standing harmonics supported in the waveguides which are 
sandwiched vertically between the temperature minimum in the lower solar atmosphere and the transition region. Standing harmonic 
oscillations, by means of solar magneto-seismology, may allow insight into the sub-resolution structure of photospheric MHD waveg- 
uides. 

Key words. Sun: atmosphere - Sun: oscillations - Sun: sunspots - Sun: photosphere 
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1. Introduction 

Over the past decades, many oscillatory phenomena have been 
observed within a wide range of magnetic waveguides in the 
solar atmosphere (B anerjee et al. 2007; Wang 2011; Asa i et all 
120121; lArreguT et al.ll2012l) . Sunspots and pores are just two of 
these many structures and they are kno wn to display sola r globa l 
oscillations, see a recent review by e.g. lPinter & Erdeivil(l201 Ih . 

The commonly studied oscillatory periods in sunspots: 3 
and 5 minutes. These oscillations are seen in intensity, line of 
sight (LOS) velocity and LOS magnetic field. The source of 
the 5-minute oscillation is a result of forcing by the 5-minute 
(/7-mode) global solar oscillation, which forms the basis of he- 
lioseismology (Thompson 2006). The 5-minute oscillations are 
typically seen in simple molecular and non-ionized metal lines, 
which form low in the umbral photosphere and are moderately 
suppressed not only in the penumbra, but also in the chromo - 
spheric atmosphere above the umbra ("Bogdan & Judge! l2006h . 
The cause of the 3-minute oscillations is still unknown but there 
are two main streams of theories: they could either be stand- 
ing acoustic w aves which are hnked to the resonant modes 
of the sunspot (iBogdan & Judge! l2006h or they could be low- 
P slow magneto-acoustic-gravity waves guided along the am- 
bient magnetic field (Bogdan 2000). The 3-minute oscillations 
are seen in plasma elements that form higher up, in the low 
chromosphere, and these are also moderately suppressed in the 



penumbra ("Christ opoulou et al.1 120001; iBogdan & Judge! 120061; 
Ikobanov et al. 200^^ 

Magnetohydrodynamic (MHD) theory, when applied to a 
cylindrical magnetic flux tube, reveals that a variety of waves 
can be supported, four of which are often reported in various 
structures in the solar atmosphere. Longit udinal (slow sausage ) 
(Ide Moortell l2009; Wang 2011), fast kink ( iAndries et al.ll2009!) . 
fast sausage ( McAteer e t al.ii2003 ) and Alfven (torsional) waves 
(iJess et al.l200 9). each of which affects the flux tube in a specific 
way. The sausage modes are of interest here; the sausage mode 
is a compressible, symmetric perturbation around the axis of a 
flux tube which causes density perturbation that can be identi- 
fied in intensity images dEdwin & Roberts!ll983h . Furthermore, 
due to the fact that the wave will either compress or expand the 
flux tube, the magnetic field will also show signs of oscillations. 
This mode may come in two forms in terms of phase speed clas- 
sification: a slow mode (often also called the longitudinal mode) 
which generally has a phase speed close to the characteristic tube 
speed and a fast mode, which has a phase speed close to the ex- 
ternal sound sp eed, assuming a region that has a plasma-p > 1 
(lErdelvill200i . The main difference between the two modes 
is the phase relationship between appropriate MHD quantities 
which allows them to be identified. In this case, the fast sausage 
mode has an out-of -phase relationship between the area and in- 
tensity, while the slow sausage mode has an in-phase relation- 
ship. The technique that was applied to obtain these phase re- 
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lationships are covered by e.g. lGoedbloed & PoedtsI (l2004l) and 
iFuiimura & Tsunetal (l2009l) . 

Sausage m o des h ave been observed in solar pores; 
iDorotovic et alJ (l2008l) observed a pore for 11 hours and re- 
ported periodicities in the range of 20-70 minutes. These os- 
cillations were consequ ently interpret e d as hn ear low-frequency 
slow sausage waves. Morton et al.l (1201 ih used the Rapid 
Oscillations in the Solar Atmosphere (ROSA) instrument to also 
identify linear sausage oscillations in a magnetic pore. However, 
determining whether the oscillations were slow or fast proved to 
be difficult. 

The source and driving mechanism(s) of these MHD sausage 
modes have been very difficult to identify. Numerical sim- 
ulations of a flux tube rooted in the photosphere are buf- 
feted by a wide range of coherent sub-photospheric drivers. 
These drivers can either be horizontal or vertical, single, paired 
or a power spectrum , wi t h varying phase differenc es (see 
e.g.,Malins & Erdelvi.2"007l: iKhomenko et al.ll2008; .Fedun et alj 
I2011allbh . One example of a horizontal driver representing solar 
global oscillations is the absorption of the global /7-mode by the 
magnetic field of the sunspot. Mathew (2008) studied this effect 
and found a structured ring-like absorption pattern in Doppler 
power close to the umbral-penumbral boundary. This effect was 
largest where the transverse magnetic field was at its greatest and 
this region allows fast waves to be converted into slow magneto- 
acoustic waves, which are a potential source of MHD waves in 
sunspots and other similar magnetic structures. 

We report here, the observation of both slow and potentially, 
fast sausage MHD waves in the lower solar atmosphere on three 
different occasions and in magnetic waveguides of various types. 
In section|2] we describe the data collection and the data process- 
ing method. In section |3] we describe the results obtained from 
the three different data series and discuss the findings. Section 
|4] details the underlying idea of identifying these oscillations as 
standing harmonics. Finally, in section|5] we conclude. 

2. Data collection and Method of Analysis 

Three time series of images with high angular resolution have 
been chosen here, in order to demonstrate the identification of 
MHD sausage waves. The images were taken in the G band 
(4305 A), which samples the low photosphere, and were ac- 
quired using: 

1 . The Swedish Vacuum Solar Te lescope (SVST) situat ed in 
La Palma on the Canary Islands. f Scharmer et al.l (Il985h pro- 
vides a detailed description of the features of the SVST. The 
images were taken on the 7th July 1999. The sunspot is in 
the Active Region NOAA 8620. The observing duration is 
133 minutes with a cadence time of 25 seconds. The field of 
view co vered an area of 33, 600 km by 54, 600 km ( 1 pixel w 
60 kml. lBonet et al.l (l2005l) gives a detailed analysis of this 
sunspot. A context image is the left image of Fig.[T] 

2. The Dutch Open Telescope (DOT), is also situated in La 
Palma on the Canary Islands. Two series of imaging data se- 
quences were taken using this telescope. A detailed guide of 
the features of the DOT is provided by Rutten et al. (2004) 
. The first series of data was taken on the 13th July 2005; 
the sunspot is in the Active Region NOAA 10789. The re- 
gion has been slowly decaying and this sunspot is leading 
the small group. The observing period is 165 minutes and 
has a cadence time of 30 seconds. The second set of data, 
taken on the 15th October 2008, is of one large pore with a 
light bridge which is about 15 pixels (750 km) wide in the 




Fig. 1: An overview of the magnetic waveguides observed for 
this analysis, (left) The 1999 sunspot observed with the SVST 
with an average umbral area of 19,650 pixels (50 Mkm^) . (mid- 
dle) The 2005 sunspot observed with the DOT with an average 
umbral area of 12,943 pixels (32 Mkm^). (right) The 2008 pore 
observed with the DOT with an average area of 10971 pixels 
(27 Mkm ^), the light bridge that separates the pore can be seen. 
Furthermore, these structures were seen near the disk centre, so 
there is little to no line of sight effects. 

Active Region NOAA 11005. The duration of the observ- 
ing run is 66 minutes and has a cadence time of 20 seconds. 
Both DOT image sequences covered an area of 50,000 km 
by 45,000km, where the maximum spatial resolution is 0.2" 
(« 140 km). Typical context images are the middle and right 
panels of Fig. [T] 

In order to obtain information relating to the cross-sectional 
area of these waveguides, a strict and consistent definition of the 
area is required - each pixel has a value of less than 3o of the me- 
dian background intensity. The background is defined as an area 
of the image where there is no formed magnetic structure. This 
may appear to be an arbitrary definition, however, a histogram 
of the background intensity reveals a Gaussian distribution and 
when adding the area around and including the waveguide, there 
is significant peak on the lower end of the Gaussian distribu- 
tion curve around 3a or higher. Thus, we have a 99% confidence 
that the area is of the structure and not of the background. The 
mean intensity value was determined by summing over the in- 
tensity of each pixel found in the waveguide and dividing it by 
the total estimated area. These waveguides are not static struc- 
tures, they slowly changed in size during the observing period. 
This slowly varying size change had to be removed in order for 
it not to mask any weaker oscillation signatures. The detrending 
was accomplished by a non-linear regression fit and the consis- 
tency of the results was compared to subtracting the residue from 
an Empirical Mode Decomposition (EMD) analysis (explained 
below). The residue is the data that remains after the EMD pro- 
cedure has extracted as many signals as possible and it provides 
a very good approximation of the background trend. 

The resulting reduced data series were then analysed with a 
wavelet tool in order to extract any periods of oscillation present 
within the data. The algorithm used is an adapted version of th e 
IDL wavelet routine developed by iTorrence & Compel (Il998l) . 
The standard Morlet wavelet, which is a plane sine wave with an 
amplitude modulated by a Gaussian function, was chosen due 
to its suitable frequency resolution. Furthermore, both white and 
red noise backgrounds were employed in order to verify the con- 
sistency of the results. The cross-hatched area marks the cone of 
influence (COI), where edge effects due to the wavelet structure 
affect the wavelet transform and anything inside the COI is dis- 
carded. The contours show the confidence level of 95%. Further 
to this, the data representing the size and intensity has also been 
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Fig. 2: Evolution of the area of the 1999 sunspot (upper panel); 
the wavelet power spectrum for a white noise background, the 
cone of influence is marked as a cross-hatched area where edge 
effects become imporant and the contour lines show the 95% 
confidence level (lower left panel). Global (integrated in time) 
wavelet power spectrum, where the dashed line shows the 95% 
confidence limit (lower right panel). 



analysed using EMD, which decomposes the time series into a 
finite number of Intrinsic Mode Functions (IMFs). These are es- 
sentially narrowband filtered time series, with each IMF contain- 
ing one or two periods that exist in the original data series. For 
more information on the features and applicability of the EMD 
method see e.g. Terradas et al. (2004). The EMD technique was 
first proposed by Huang et al.l (119981) and offers certain benefits 
over more traditional methods of analysis, such as wavelets or 
Fourier transforms. Generally, the next step after EMD analy- 
sis is to construct a Hilbert power spectrum which has a better 
time and spatial resolution than either wavelet or EFT routines. 
However, this has not been carried out due to a lack of a robust 
code base at this time and will be addressed in future work. At 
this stage, we rely on wavelet and EMD analyses, as customary 
in solar physics. 



3. Results and Discussion 

3.1. Sunspot, 7 July 1999, AR 8620 

Fig.|2]shows the wavelet analysis of the 1999 sunspot area data. 
There are three confidently identified periods that exist with 95% 
certainty; 5, 14 and 32 minutes. The 32-minute period is found 
over a wide range of the time series, with some of its power in- 
side the COL However, most is confidently outside the COL The 
14-minute period is highly localised at around 50 to 120 min- 
utes of the data series and it appears to be an off shoot of the 32- 
minute period. The 14-minute period starts at 17 minutes and it 
slowly decays and stabilizes at 14 minutes. The 5-minute period 
is relativity weak compared to the other two periods and there 
are regions of strong power where the period is either greater or 



lower than 5 minutes and this could be indicative of the 3-minute 
oscillation. 

The intensity wavelet shows three distant periods of oscilla- 
tions: 6, 24 and 33 minutes. The 6- and 33-minute periods both 
correspond to two of the area wavelet oscillations, however, the 
24-minute oscillation is not seen in the area wavelet. Some effect 
must therefore be causing the disparity in oscillation periods, 
perhaps the opacity effect (see Fujimura & Tsuneta 2009). It is 
safe to say that these oscillations are caused by sausage waves. 
The reason is that in standard MHD theory, the sausage wave 
is the only MHD wave capable of ch anging the area of th e flux 
tube that is observed on disk (see e.g. Cooper et al.ll2003"allbl) . 

Without the ability to directly compare the area to intensity, 
great caution needs to be exercised to determine with confidence 
whether the perturbations are fast or slow. A wavelet phase dia- 
gram reveals regions (where the wavelet coherence is high and 
the period is < 20 minutes) to be either out-of -phase or in-phase 
but a clear image of constant phase difference does not appear. 
This might be due to mode conversion occurring in the sunspot, 
since the G-band samples a region where plasma-P is « 1 in a 
magnetic structure (Garv 2001). When the period is > 20 min- 
utes, the only area of high coherence is located around 32 min- 
utes and found to be out-of-phase, which hints that there is a fast 
sausage wave. However, only two full wave periods are outside 
the COI, which is due to the total length of the data series. 

Fig. |3] shows the computed IMFs for the 1999 sunspot data 
set. The IMFs show the periods of oscillations identified using 
the EMD routine. Six IMFs are shown, two were neglected due 
to uncertainties and the additional residue is ignored. In gen- 
eral, the higher order IMFs tend to show longer periods and as 
such contain fewer wave periods, which makes phase identifica- 
tion less reliable. In this case, only two IMFs coincide (c3 and 
ce) with the wavelet period that shows both area and intensity 
perturbations. IMF C3 shows a period of 5 minutes and contains 
several regions of three or more wave periods that are either in- 
or out-of-phase behaviour which agrees with the cross-wavelet 
phase analysis. IMF ce has clear out-of-phase behaviour between 
the area and intensity for first three wave periods, which is po- 
tentially indicative of a fast MHD sausage wave. 

It was possible to approximately separate the penumbra from 
the umbra and investigate its area for oscillations. However, the 
penumbra is a highly dynamic object and this makes the area 
estimation reasonably uncertain. There seem to be four periods 
that exist at 95 % certainty: 7, 14, 25 and 44 minutes. The longest 
period happens to have the majority of its power inside the COI, 
so it has been discounted. The 26-minute period is mostly con- 
centrated at the start of the time series but as time evolves for- 
ward, the power starts to decay and the period drops to 20 min- 
utes. This change in period is most likely caused by the decay 
of the sunspot over the observational period and the fact that 
the magnetic flux slowly decays along with it. The two shorter 
periods (7 and 14 minutes) closely correspond to the 5- and 14- 
minute oscillations in the umbra; they could be a continuation 
of these umbral periods that became up-shifted as they enter the 
less compact structure of the penumbra. The EMD analysis of 
the penumbral oscillations reveals that area oscillates in-phase 
with the intensity. No signal of fast sausage MHD wave is seen; 
this indicates in this case that mode conversion has not occurred 
in this region. 

3.2. Sunspot, 13 July 2005, AR 10789 

Fig. |4] shows the wavelet analysis of the 2005 sunspot area in 
AR 10789. There are six periods that exist at 95 % confidence 
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Fig. 3: The IMFs of the evolution of the area (red) and intensity (black) for the 1999 sunspot, over-plotted to aid comparison. 
Generally after the 6th IMF, higher IMFs lack a sufficient number of wave periods, which makes it difficult, and less reliable, to 
obtain an accurate period. 
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Fig. 5: Same as Fig.[3]but for the sunspot in Active Region 10789 in 2005. 
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Fig. 4: Same as Fig.[2]but for the sunspot in Active Region 10789 
in 2005. 



Fig. 6: Same as Fig. |2]but for the pore in Active Region 1 1005 
in 2008. 



level: 4, 5.5, 7.5, 10.5, 17.5 and 65 minutes. The longest period 
has its wavelet power reach a maximum value at 65 minutes but 
the power is spread from 50 to 80 minutes and, in this case, it is 
inside the COl and so it has been neglected for further analysis. 
A comparison with the corresponding intensity wavelet reveals 
that the intensity shows 7.5, 10.5, 21 and 65-minute oscillations. 
The cross-wavelet phase indicates that these oscillations are in- 
phase. There appears to only be very small and infrequent re- 
gions of out-of-phase behaviour that occur for 1 to 2 wave peri- 
ods at around periods of 2 to 3 minutes. This is most likely an 
artefact of the wavelet transform at very low periods. 

Fig. |5] shows the IMFs for the area and the intensity of the 
sunspot data in AR 10789. Only C3 and C4 show oscillations that 
occur in the area and the intensity time series. The region of 
interest is within the time interval of 90-130 minutes, where the 
wavelet has these oscillations. Both IMFs show clear in-phase 
behaviour in this time interval and both have one wave period 
out-of-phase at around 110 minutes. The overall phase relation 
between the area and intensity indicates the presence of slow 
sausage waves. 

3.3. Pore, 1 5 October 2008 

Fig.|6]shows the wavelet analysis of the pore with a light bridge. 
There are three periods that exist at 95 % confidence level: 4.5, 
8.5 and 14.5 minutes. The majority of the power of the period 
of 14-15 minutes is inside the COI and so this period has been 
discarded. The other two periods: 4.5 and 8.5 minutes are seen in 
both area and intensity data when the wavelet analyses are cross- 
correlated. The power for these two periods is concentrated in 
the time interval of 20-60 minutes. The cross-wavelet analysis 
shows that the overlapping time span is somewhat smaller, at 
about 30-50 minutes. 

Fig. I2] shows the IMFs for the area with intensity over- 
plotted. In this case, C3 indicates a period of 4.5 minutes and 
C5 has a characteristic period of 8.5 minutes; this applies to both 



the area and intensity IMFs. C3 reveals that the phase relation 
is out-of-phase for about three wave periods, with a period of 
in-phase behaviour only at either end of this group for the time 
interval 20-40 minutes. This behaviour occurs twice more in this 
IMF. The comparison of C5 also reveals that the phase relation 
is out-of-phase for about 2 wave periods and shows a small re- 
gion of in-phase behaviour at the end of this time series. Once 
again, this behaviour is a potential indicator of the presence of 
fast sausage MHD waves in this pore. 

Let us make a few general comments apply to all the results. 
The easiest way to confirm the linearity of waves is to compare 
the amplitude of the oscillations to the characteristic scale of the 
structure. In all three cases studied here, the oscillation ampli- 
tudes are around 10% or less of the total area, which indicates 
that these oscillations are linear. Furthermore, the amplitude of 
the oscillation in these three cases is roughly the same. So the 
amplitude scales with the size of the structure. 

4. Standing Harmonics 

Basic MHD theory interpretation allows sunspots and pores to 
be described as vertical cylindrical flux tubes, with the base 
bounded in the photosphere and the top bounded at the transi- 
tion region due to the sharp gradients in the plasma properties at 
these locations dLuna-Cardozo et al]|2012h . Taking this further, 
an ideal flux tube is assumed here. The plasma density and mag- 
netic field are homogeneous within the flux tube. This means that 
the standing harmonics of such flux tubes are the MHD equiv- 
alent to those of the harmonics in an open-ended compressible 
air pipe, where the ratio of the harmonic periods is given by. 
Pi /Pi ^2, Pi /Pi = 3 and so forth. 

Let us now summarise the observed findings. Table [T] con- 
tains the periods of oscillations found in all three magnetic 
waveguides. 

For the 1999 sunspot, there are 3 periods found. The second 
period at 14 minutes gives a period ratio (P1/P2) of 2.3 ± 0.3, 
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Fig. 7: Same as Fig.[3]but for the pore in Active Region 11005 in 2008. 



Data Set 


Period (Mins) 


Ratio (Pi /Pi) 




Pi 


-32 ±3 






Sunspot 1999 


Pi 


- 14± 1 


2.3 


±0.3 




Pi 


-5 ±0.5 


6.4 


±0.9 




Pi - 


17.5 ± 1.5 






Sunspot 2005 


Pi- 


10.5 ± 0.5 


1.7 


±0.2 




Pi- 


7.5 ±0.5 


2.5 


±0.3 




Pa- 


5.5 ±0.5 


3.2 


±0.4 




Ps 


- 4 ± 0.5 


4.4 


±0.7 


Pore 2008 


Pi - 
Pi- 


8.5 ±0.5 
4.5 ±0.5 


1.8 


±0.2 



Table 1 : The periods of oscillations that are found in the area of 
the waveguides that exist at 95% confidence level and are out- 
side the COL There is another period for the 2008 pore which 
is discussed in the text but not included in the table due to the 
criteria laid out. 



which is slightly higher than the expected value of a uniform 
waveguide with a canonical value of 2 but it is within the error 
range. The other period is more controversial to fit into this im- 
age. The 5 minute period has a ratio that is maybe too high to 
incorporate into a harmonic standpoint and begs the question: 
where are the other harmonics in-between? It is more plausible 
that this 5 minute period is just a result from the global p-mode 
buffeting the magnetic waveguide. 

For the 2005 sunspot in AR 10789, there is a clearer picture 
of potential harmonics. The first period is at 17.5 minutes and 
the second period is at 10.5 minutes, which gives a ratio of 1.7 
± 0.2, and the third period at 7.5 minutes gives a ratio of 2.5 
± 0.3. The period ratio is modified downwards in a consistent 



manner, by 0.2 or 0.3 as the harmonic number increases. These 
ratios are strong evidence for standing waves in this magnetic 
waveguide. The last two periods (5.5 and 4 minutes) are again 
too close to the global p-mode oscillation period. However, the 
5.5 minute period still fits the pattern for the decreasing period 
ratio whereas the last period breaks this pattern. 

For the 2008 pore of AR 1 1005, the picture is more muddled 
due to the short time series available. The 15 minute period has 
been ignored due to having the majority of its power inside of 
the COl but is more likely to be the first harmonic than the 8.5 
minute period. Assuming that the 8.5 minute period is the first 
harmonic, the ratio is 1.8 for the 4.5 minute period, which hap- 
pens to be very similar to the ratio of the 2005 sunspot's first 
period ratio. If we ignore the issue with the 14.5 minute period 
for the time being, it alters the period ratios to 1 .7 and 3.2 for 8.5 
and 4.5 minutes respectively. This result makes more sense as it 
is more likely that any observation of periods around 5 minutes 
is a direct result of the global p-mode oscillation and this brings 
the second period ratio into line with the 2005 sunspot period 
ratio. 

The main conclusion to take away from this data analysis so 
far is that the simple homogeneous flux tube model cannot fully 
account for these ratios. However, this simple model seems to 
be robust enough to give a good first insight. The most likely 
reasons for deviation from the canonical period ratio value are 
firstly that sunspots and pores (just like most lower atmospheric 
magnetic structures) expan d wit h height, causing magnetic strat- 
ification (IVerth & Erdelvill2008tlLuna-Cardozo et al.ll2012l) . and 
secondly, that the Sun's gravity causes density stratification 
(Andries et al. 2009). These two effects will either increase or 
decrease the period ratio of the harmo nics depending on the cho- 
sen density or magnetic profile (see iLuna-Cardozo et al ] l2012i 
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for a detailed analysis in the context of slow sausage oscilla- 
tions). In addition, these magnetic structures are rarely purely 
cyhndrical, they ca n be elli ptical (or arbitrary) in shape (see 
iRuderman & Erde lvi 2009; Morton & Erdelvil2009l) and in most 
cases are non-axially symmetric. Also, in some cases the flux 
tube is more suitably described as open-ended at the transition 
region, which would remove the even harmonics. 

5. Conclusions 

In this paper we have investigated several magnetic waveguides, 
with the objective of detecting MHD sausage waves and deter- 
mining whether they are slow or fast, propagating or standing. 
Based on the results presented here, we have confidently inter- 
preted the observed periodic changes in the area cross section 
of flux tubes, which are manifested as a pore and two sunspot 
waveguide structures, as proof of the existence of linear slow 
and fast sausage MHD oscillations. Using wavelet analysis, we 
found standing waves in the photosphere with periods ranging 
from 4 to 32 minutes. Employing complementary EMD anal- 
ysis has allowed the MHD modes detected to be identified as 
a combination of fast sausage and slow sausage modes, due to 
the phase behaviour of the area and intensity. It is very likely 
that these oscillations are standing harmonics supported in a flux 
tube. The period ratio {P\/Pi=23) of these oscillations indicates 
strongly that they are part of a group of standing harmonics in 
a flux tube that is non-homogeneous and is bound by the photo- 
sphere and the transition region. Furthermore, there is possible 
indirect evidence of mode conversion occurring in one of these 
magnetic waveguides. 
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